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This  study  investigates  the  dynamic  impacts  of  GDP  growth,  energy  consumption  and  population  growth 
on  C02  emissions  using  econometric  approaches  for  Malaysia.  Empirical  results  from  ARDL  bounds 
testing  approach  show  that  over  the  period  of  1970-1980,  per  capita  C02  emissions  decreased  with 
increasing  per  capita  GDP  (economic  growth);  however  from  1980  to  2009,  per  capita  C02  emissions 
increased  sharply  with  a  further  increase  of  per  capita  GDP.  This  is  also  supported  by  the  dynamic 
ordinary  least  squared  (DOLS)  and  the  Sasabuchi-Lind-Mehlum  U  (SLM  U  test)  tests.  Consequently,  the 
hypothesis  of  the  EKC  is  not  valid  in  Malaysia  during  the  study  period.  The  results  also  demonstrate  that 
both  per  capita  energy  consumption  and  per  capita  GDP  has  a  long  term  positive  impacts  with  per  capita 
carbon  emissions,  but  population  growth  rate  has  no  significant  impacts  on  per  capita  C02  emission. 
However,  the  study  suggests  that  in  the  long  run,  economic  growth  may  have  an  adverse  effect  on  the 
C02  emissions  in  Malaysia.  Thus,  significant  transformation  of  low  carbon  technologies  such  as 
renewable  energy  and  energy  efficiency  could  contribute  to  reduce  the  emissions  and  sustain  the  long 
run  economic  growth. 
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1.  Introduction 

The  global  concern  for  environmental  sustainability  is  apparent 
in  an  increasing  public  mandate,  and  the  development  strategy 
largely  depends  on  whether  the  enduring  economic  growth  causes 
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environmental  degradation  or  whether  such  growth  is  sufficient  to 
compensate  for  the  environmental  cost  of  production  or  develop¬ 
ment  process.  Nevertheless,  exhaustible  natural  resources  serve  as 
inputs  into  the  production  or  development  process.  If  the  functional 
relationship  between  natural  resources  and  modern  development 
processes  cannot  be  avoided,  then  harm  to  the  environment  is 
inevitable  [1].  The  prevalence  of  such  problems  is  higher  in 
countries  such  as  Malaysia,  where  economic  growth,  energy  security 
and  environmental  sustainability  are  simultaneously  important. 


http://dx.doi. org/10.1016/j.rser.2014.07.205 
1364-0321/©  2014  Elsevier  Ltd.  All  rights  reserved. 
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In  Malaysia,  the  government  has  set  a  voluntary  target  to  reduce 
40%  carbon  emission  intensity  [2]  as  well  as  a  goal  of  achieving 
high  income  nation  by  2020  [3],  Thus,  a  significant  concern  is 
whether  the  goals  of  higher  economic  growth  and  improved 
environmental  quality  (emission  reduction)  are  mutually  exclu¬ 
sive.  For  instance,  Jorgenson  and  Wilcoxen  [4]  observed  that  slow 
economic  growth  coincided  with  the  advent  of  environmental 
regulations  when  petroleum  prices  increased  in  the  United  States 
of  America  (USA).  However,  the  Malaysian  economy  subsidises 
petroleum  prices  [5]  which  alternatively  encourages  additional 
energy  consumption.  Such  subsidised  prices  might  be  conducive 
for  accelerating  economic  growth.  Moreover,  from  1970  to  1980, 
Malaysia's  economic  dominance  shifted  from  the  primary  sector 
(mainly  agricultural)  to  the  secondary  sector  (industrial)  [6]; 
consequently,  the  energy  consumption  rose  substantially.  Cur¬ 
rently,  the  economic  development  in  Malaysia  is  moving  towards 
the  services  sector,  which  is  a  highly  energy  intensified  sector. 
Therefore,  achieving  a  reduction  of  40%  emission  intensity  might 
be  distorted,  unless  the  economy  substantially  shifts  to  the  use  of 
low  carbon  technology  (renewable  energy  and  energy  efficient)  to 
produce  the  goods  and  services.  Although  the  theory  of  the 
environmental  Kuznets  curve  (EKC)  refers  to  an  early  stage  of 
economic  development,  whereby  a  small  portion  of  the  excess 
income  achieved  by  industrialisation  is  positively  associated  with 
environmental  problems  [1],  Therefore,  this  study  aims  to  inves¬ 
tigate  the  dynamic  impacts  of  GDP  growth,  energy  consumption 
and  population  growth  on  the  C02  emissions  in  Malaysia  by 
applying  the  ARDL  bounds  testing  approach  and  dynamic  OLS. 
Moreover,  the  Sasabuchi-Lind-Mehlum  (SLM)  U  Test  and  the 
Quadratic  and  Polynomial  Confidence  Interval  diagram  are  also 
performed  to  check  the  robustness. 

2.  Review  of  literature 

A  number  of  empirical  studies  investigated  the  relationship 
between  C02  emissions  and  economic  growth  for  a  group  of 
developed  and  or  developing  countries  such  as,  BRIC  countries  [7]; 
MENA  countries  [8];  OECD  countries  [9];  industrialised  countries 
[10];  and  small-economy  countries  [11],  Most  of  these  studies 
especially  on  the  EKC  provide  a  simple  path  of  economic  growth 
and  environmental  problems  such  as  emissions  or  pollutions.  Some 
studies  also  focused  this  relationship  on  a  single  country  for  instance, 
Bangladesh  [12];  China  [13-19];  Turkey  [20];  South  Korea  [21,22]; 
India  [23];  Nigeria  [24];  and  the  USA  [25,26],  The  relationship 
between  C02  emission  and  economic  growth  may  appear  a  linear, 
U  shape,  inverted  U  shape  or  any  other  shapes  [27], 

However,  some  studies  confirmed  the  EKC  hypothesis  with 
an  inverted  U-shaped  relationship  [28-34,21]  which  indicates  an 


additional  increase  of  economic  growth  improves  the  environmental 
quality  i.e.  emission  reduction.  In  contrast,  many  studies  found  non¬ 
existence  of  EKC,  indicates  increasing  economic  growth  leads  a 
greater  environmental  degradation  (e.g.  [35-40]).  However,  Holtz- 
Eakin  and  Selden  [41]  found  a  monotonic  rising  curve,  whereas 
others  found  a  non-monotonic  relationship  [42-45],  Furthermore, 
Shafik  [44]  and  Grossman  and  Krueger  [1]  found  an  N-shaped  cun/e 
while  Ozturk  and  Acaravci  [46]  found  no  causal  relationship  between 
C02  emission  and  GDP  per  capita.  In  addition,  a  few  studies  found  a 
unidirectional  causality  from  energy  consumption  to  aggregate  out¬ 
put  or  income  [47-49]  whereas  other  studies  found  bi-directional 
causality  between  energy  consumption  and  aggregate  output  [50— 
55],  Nevertheless,  each  of  these  studies  has  its  own  merits  for  better 
understanding  of  the  paradox  of  economic  growth  and  environ¬ 
mental  balance. 

In  Malaysia,  there  is  a  veiy  limited  study  on  the  dynamic 
impacts  between  economic  growth,  C02  emissions  and  other 
variables.  Table  1  summarises  their  findings,  analytical  techniques 
and  limitations.  A  few  studies  found  a  positive  and  bidirectional 
relationship  whereas  Saboori  et  al.  [56]  found  an  existence  of  EKC 
for  Malaysia  from  1980  to  2009.  Unlike  these  relationships,  a  few 
studies  also  demonstrated  the  evidences  of  C02  emission  reduc¬ 
tion  with  policies  such  as  public  transport  and  waste  to  green 
energy  [57-59], 

Based  on  the  above  literatures,  the  empirical  results  demon¬ 
strated  a  varied  relationship.  This  might  be  due  to  the  selection  of 
different  samples  (a  country/few  country/a  region);  time  periods; 
variables  and  analytical  techniques.  However,  this  study  would  be 
the  first  attempt  to  investigate  the  dynamic  relationships  between 
C02  emissions  and  other  variables  such  as,  energy  consumption, 
GDP  per  capita  and  population  growth  in  Malaysia.  The  energy 
consumption  and  population  growth  includes  as  explanatory 
variables  because  omitted  variables  may  often  produce  misleading 
results  from  the  OLS  estimator  [60]  which  also  helps  to  fill  the 
research  gaps. 
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Fig.  1.  Logarithmic  trend  of  per  capita  C02  emissions,  GDP  and  energy  consump¬ 
tion;  and  population  growth  rate. 


Table  1 

Summary  of  studies  on  economic  growth,  C02  emissions  and  other  variables  in  Malaysia. 


Relationships  and  sources 


Techniques 


Limitations3 


Positive  relationship  between  economic  growth  and  energy 
consumption,  and  emissions;  Ang  [61] 

Predicted  C02  emissions  and  energy  consumption  should  be  tripled 
to  maintain  the  GDP  at  4.6%  by  2030;  Gan  and  Li  [62] 

Bidirectional  relationship  between  C02  emissions  and  economic 
growth;  Azlina  and  Mustopha  [63] 

C02  emissions  decrease  with  an  increase  of  GDP  per  capita,  which 
indicates  the  existence  of  EKC;GDP  per  capita  and  its  quadratic 
form  as  explanatory  variables;  Saboori  et  al.  [56] 


Johansen  cointegration; 
VECM-based  Granger 
causality  test 

Ordinary  least  squares  (OLS) 
estimation 

Johansen  cointegration  test; 
VECM-based  Granger 
causality  test 
ARDL  bounds  testing  to 
Cointegration;  VECM-based 
Granger  causality  test 


Did  not  consider  the  quadratic  form  of  GDP  per  capita  to  examine 
the  existence  of  EKC. 

No  unit-root  test,  although  macroeconomic  variables  are  likely  to 
be  characterised  by  unit  root  processes  [64]. 

Presence  of  mixed-order  integration.  Johansen  technique  is  highly 
sensitive  to  the  selection  of  optimal  number  of  lags  [65]. 

Ignore  other  variables,  such  as  energy  consumption  and  population. 
Inclusion  of  the  quadratic  form  of  the  GDP  per  capita  fulfils  the 
necessary  conditions  for  the  existence  of  an  inverted  U-shaped 
relationship  but  not  sufficient  conditions  [66]. 


Author's  view. 
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3.  Outlook  of  Malaysia's  economy,  energy  consumption  and 
emissions  to  2030 

Gan  and  Li  [62]  estimated  the  expected  GDP  growth  of  Malaysia 
averaged  4.6%  between  2004  and  2030,  and  the  actual  GDP  would 
reach  US$341.6  billion  by  2030  compared  to  US$107  billion  in 
2004.  A  structural  change  may  occur  in  Malaysia  for  instance,  the 
output  of  agriculture  sector  might  drop  by  3%,  whereas  the 
industrial  and  service  sectors  might  marginally  increase  over  the 
period  of  2004-2030.  In  addition,  it  is  anticipated  that  the  per 
capita  real  GDP  would  reach  US$8690  in  2030  from  US$  4296  in 
2004.  However,  Gan  and  Li  [62]  also  projected  that  energy 
consumption  would  also  increase  by  4.3%  in  2030  to  keep  pace 
with  economic  growth.  Although  the  growth  rate  of  energy 
consumption  is  a  bit  lower  than  the  GDP  growth  rate,  the  absolute 
amount  of  energy  consumption  would  be  3  times  greater. 

Fig.  1  illustrates  the  logarithmic  form  of  C02  emissions  per 
capita,  GDP  per  capita,  energy  consumption  per  capita,  and 
population  growth  rate.  It  clearly  demonstrates  that  the  C02 
emissions  per  capita,  GDP  per  capita,  and  energy  consumption 
per  capita  are  constantly  increasing,  whereas  the  population 
growth  rate  is  decreasing.  The  rising  intensity  of  C02  emissions 
is  of  a  great  concern  due  to  the  Malaysia's  rapid  economic  growth. 

In  Malaysia,  power  sector  (electricity)  is  the  major  source  of 
energy  consumption  and  C02  emissions.  This  sector  accounted  for 
approximately  36%  of  the  total  global  C02  emissions  from  oil  in 
2004  [67].  Among  all  the  greenhouse  gases,  C02  emissions  are 
responsible  for  more  than  60%  of  the  greenhouse  gas  effects, 
which  cause  to  global  warming  and  climate  change.  Therefore, 
mitigation  or  reduction  of  C02  emissions  became  a  global  concern 
to  secure  sustainable  energy  as  well  as  to  minimise  the  adverse 
effects  of  climate  change. 


4.  Data  and  analysis 

The  study  used  annual  time  series  data  from  1970  to  2009  for 
Malaysia,  which  were  obtained  from  the  world  development 
indicator  (WDI)  dataset.  The  variables  of  interest  included  C02 
(carbon  dioxide)  emissions  per  capita  (metric  tons  per  capita)  as  a 
dependent  variable  and  GDP  per  capita  (GDPC),  energy  consump¬ 
tion/use  (kg  of  oil  equivalent  per  capita)  (EUC)  and  population 
growth  (POPG)  as  explanatory  variables.  The  rational  for  specifying 
the  employed  model  and  variables  chosen  discussed  as  follows. 

The  empirical  model  of  this  study  derives  by  the  following 
standard  Cobb-Doglus  production  function  with  constant  returns, 
the  aggregate  output  function  can  be  displayed  at  time  t  as  follows: 

Yt  =  F(Kt,ALt) 

where  Yt  is  GDP,  Kt  is  capital  and  ALt  is  the  effective  labour. 

Since  it  is  widely  assumed  that  C02  emission  is  discharged  by 
economic  activities,  therefore  C02  emission  function  can  be 
written  as: 

CO  2(t)  =  p(F(Y(t)) 

where,  v  represents  a  certain  rate  of  C02  emission  from  the 
production  function. 

In  the  production  function,  all  forms  of  capital  are  not  respon¬ 
sible  for  C02  emission  while  burning  of  energy,  e.g.  oil,  coal,  gas 
and  electricity  are  mainly  responsible  for  discharging  C02  emis¬ 
sion.  Hence,  total  capital  can  be  the  composition  of  emitting 
capital  ( I(e )  and  non-emitting  capital  ( Kn )  as  shown  below 

K  =  I<e+Kn 


Therefore,  C02  emission  function  can  be  rewritten  as: 

C02(t)  =  4>Ke(Y) 

Since,  Y  is  the  function  of  labour  and  capital;  hence,  population 
growth  can  be  taken  as  a  proxy. 

Moreover,  the  concentration  of  emissions  is  increasing  as  a 
consequence  of  various  human  activities  where  population  growth 
has  been  a  core  factor  to  explain  C02  emission  dynamics  [68], 
Thus,  the  model  can  be  written  as: 

In  COt2  =  /?o +/h  GDPt  +/?2Kpt  +PoPCt  +  et 

Since,  Grossman  and  Krueger  (1995)  have  shown  a  non-linear 
relation  between  GDP  and  C02  emission  that  can  be  reflected  by 
the  following  equation: 

In  COt2  =/?„+/?!  In  GDPt+/72  In  GDPt  +/?3  In  I<pt+fi4PGt  +  et  (1) 

Hence,  this  study  employed  the  autoregressive  distributed  lag 
(ARDL)  technique  developed  by  Pesaran  et  al.  [69]  to  assess  the 
long-run  relationship  and  the  changes  of  the  key  variables.  The 
key  characteristics  of  this  technique  include  (i)  a  cointegration 
relationship  estimated  using  OLS  and  estimation  performed  after 
selecting  the  respective  lag  order  appropriate  for  the  model;  (ii)  an 
applied  technique  that  remains  statistically  significant,  regardless 
of  the  nature  of  the  variable  integration  (i.e.,  either  1(0 )  or  /(l)  or 
mutual  co-integrated)  that  typically  explains  the  position  that  the 
unit  root  tests  may  not  be  necessary,  notwithstanding  the  Johan¬ 
sen  and  Jeslius  approaches;  and  (iii)  a  requirement  that  the  test  is 
valid  for  a  small  and  finite  data  size.  Based  on  Eq.  (1),  the  ARDL 
version  of  the  vector  error  correction  model  (VECM)  can  be 
specified  as  Eq.  (2). 

A  In  C02  —  Fo  F  F\  In  C02t  —  i  ~f-F2LGDPC  f — i  T F3  GDPGSf  _  i  -t- /^POPGf  —  l 

f  Yi  In  C02t-i+  ^  <5jAGDPGt_j 
'  =  i  j  =  i 

+  2  <p/AGDPSf_/-{-  f?mPOPGt-m+  X  ^r^UCf _r 

(=1  m  =  1  r  =  1 

(2) 

The  estimation  procedure  begins  with  Eq.  (2)  and  uses  OLS  to 
enable  the  Wald  test  or  F-test  to  determine  the  joint  significance  of 
coefficient  of  the  lagged  variables.  The  essence  of  this  procedure  is 
to  examine  the  likelihood  of  any  possible  long-run  relationship 
among  the  respective  variables.  In  this  regard,  the  null  hypothesis 
(H0)/h  =  p2  =  P3  =  Pa  =  Ih  =  0  explains  that  cointegrating  relation¬ 
ships  do  not  exist  among  the  regressors  and  the  regressand: 
therefore,  the  alternate  hypothesis  is  (Ha)  /A  ^  P2  ^  P3  ^  /?4  ¥= 
/?5  =£  0.  The  F-statistics  can  be  compared  against  the  critical  values 
of  the  upper  and  lower  bounds,  as  in  Pesaran  et  al.  [69],  Thus,  the 
whole  notion  of  this  arrangement  is  that  if  the  F-statistics  are 
higher  than  the  upper  critical  value,  then  the  null  hypothesis 
shows  that  a  cointegrating  relationship  does  not  exist.  So,  the  null 
hypothesis  is  rejected.  This  means  that  the  long-run  relationship 
exists  among  the  respective  variables,  whereas  if  the  F-statistics 
are  less  than  the  lower  critical  value,  the  null  hypothesis  is 
accepted.  Alternatively,  if  the  F-statistics  are  observed  to  be  within 
the  lower  and  upper  critical  values,  then  the  test  is  inconclusive. 
When  this  procedure  is  completed,  the  next  step  is  to  proceed 
with  the  estimation  of  the  long-run  coefficient  of  the  ARDL  model 
using  Eq.  (3). 

In  C02  =  fi0  +  2  Yi  In  C02t_,  £  <5jGDPGt_j  +  £  <p(GDPGS2t_/ 

i  =  1  j  =  0  1  =  0 

~F  ^  /7mPOPGt_m+  t//r  In  EUCf_r  +  £t  (3) 

m  =  0  r  =  0 

This  study  also  estimates  the  error  correction  model  presented  in 
Eq.  (4)  to  investigate  the  short-run  dynamics  of  the  respective 
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variables  along  with  the  short-run  adjustment  rate  towards  the 
long-run  rate. 

A  In  CO2  =  /?o  +  2  j'jA  In  C02t_i  X  <5jAGDPGt_j+  ^  <?,/AGDPGS2t-/ 

i  =  1  j  =  1  /  =  1 

4-  ^  ^/mAPOPGf  _m  +  ^  vg-A  In  EUCt_r-b  #emCf_  1  +  et 

m  =  1  r  =  1 

(4) 

Finally,  this  study  examines  the  stability  tests  of  the  coefficients  by 
employing  the  cumulative  sum  of  recursive  residuals  (CUSUM) 
and  cumulative  sum  of  squares  of  recursive  residuals  (CUSUMSQ) 
according  to  the  method  by  Pesaran  and  Pesaran  [70],  In  addition, 
this  study  also  conducts  the  Sasabuchi-Lind-Mehlum  (SLM)  U  test 
to  ensure  the  inverted  U-shaped  relationship.  It  should  be  noted 
that  Lind  and  Mehlum  [66]  created  and  altered  the  likelihood  ratio 
test  of  Sasabuchi,  which  is  referred  to  as  the  Sasabuchi-Lind- 
Mehlum  (SLM)  U  test. 

5.  Results  and  discussion 

As  mentioned  in  the  methodology  section,  prior  to  the  estima¬ 
tion  of  cointegration,  it  is  not  necessary  to  check  the  order  of 
integration  of  the  respective  variables  for  the  methodology  of 
Pesaran  et  al.  [69],  However,  this  study  conducted  the  unit  root  for 
two  reasons:  first  to  ensure  that  no  variable  surpassed  the  order  of 
integration  /( 1)  and  then  to  justify  the  appropriateness  of  applying 
the  ARDL  approach  rather  than  the  standard  cointegration 
approaches.  Therefore,  the  Dickey-Fuller  generalised  least  squared 
(DFGLS)  approach  is  employed  based  on  the  trends  and  constants 
to  examine  the  autoregressive  unit  root  by  Elliot  et  al.  [71],  The 
DFGLS  is  a  simple  modified  version  of  the  conventional  augmented 
Dickey-Fuller  (ADF)  t-test  that  de-trends  the  series  prior  to  the 
estimation  of  ADF  test  regression. 

Table  2  shows  that  the  logarithmic  form  of  C02  emissions  (per 
capita)  was  stationary  at  1(0).  However,  the  LGDPC  was  found  to 
non-stationaiy  at  /(0),  but  it  was  stationary  after  taking  first 
difference.  The  DFGLS  test  also  revealed  that  the  population 
growth  rate  was  non-stationary  at  1(0)  but  stationary  at  /(l). 
Similarly,  the  energy  use  (per  capita)  was  stationary  at  /(l)  rather 
than  at  1(0).  Therefore,  the  presence  of  such  mixed  orders  of 
integration,  as  reported  in  the  DFGLS  test,  endorsed  the  applica¬ 
tion  of  the  ARDL  approach  instead  of  conventional  econometric 
approaches.  Before  estimating  Eq.  (2),  the  optimum  lag  has  been 
identified.  Based  on  the  results  reported  in  Table  3,  the  study 
selects  the  optimum  lag  to  be  3,  according  to  the  Schwarz 
information  criterion. 

To  examine  the  cointegrating  relationship,  Wald  test  or  F-test 
for  the  joint  significance  of  the  coefficient  of  the  lagged  variables 
was  applied.  Therefore,  the  null  hypothesis  (Ho)/?]  =  f)2  = 
p3=p4  =  0  explained  the  existence  of  no  cointegration  against 
the  alternative  hypothesis  of  px  ^  p2  p3  ^  /J4  ^  0,  which  showed 
that  cointegration  existed.  Table  4  reports  the  results  of  the 


Table  2 

Results  from  Dickey-Fuller  generalised  least  squared  (DFGLS)  unit  root  test  with 
trends  and  constants. 


Log  levels  (Zt) 

Log  1st  difference  (Zt) 

/(d) 

Variable 

SIC  lag 

DFGLS  stat 

Variable 

SIC  lag 

DFGLS  stat 

lco2 

1 

-2.306** 

alco2 

0 

-8.737*** 

7(0) 

LGDPC 

0 

-1.684*** 

ALGDPC 

1 

-4.873*** 

7(0) 

POPG 

3 

-2.061 

APOPG 

2 

-4.114** 

/(l) 

EUC 

0 

-2.331 

AEUC 

0 

-6.083*’* 

/(l) 

Note:  Statistically  significant  at  1%  (***),  5%  (**)  and  10%  (*)  level. 


Table  3 

Lag  length  selection  criteria  for  cointegration. 


Lag 

Log  L 

LR 

FPE 

AIC 

SC 

HQ 

0 

-340.310 

NA 

4119.191 

19.674 

19.852 

19.736 

1 

- 194.502 

249.954 

2.493 

12.257 

13.146 

12.564 

2 

-149.216 

67.279 

0.486 

10.583 

12.183 

11.136 

3 

- 108.263 

51.486" 

0.129 

9.157 

11.468" 

9.955 

4 

-92.448 

16.266 

0.161 

9.168 

12.190 

10.211 

5 

-64.463 

22.387 

0.119" 

8.483" 

12.216 

9.772" 

LR:  sequential  modified  LR  test  statistic  (each  test  at  the  5%  level)  FPE:  final 
prediction  error,  AIC:  Akaike  information  criterion,  SC:  Schwarz  information 
criterion,  HQ.:  Hannan-Quinn  information  criterion. 

a  Indicates  lag  order  selected  by  the  criterion. 

calculated  F  statistics  in  which  each  variable  was  normalised  as 
the  dependent  variable. 

The  calculated  F  statistic  FCo2(C02ILGDPC,  LGDPCS,  POPG,  EU)= 
5.444  is  greater  than  the  upper  bound  critical  value  of  Pesaran 
et  al.  [69]  at  the  1%  significance  level  (4.15).  This  result  indicates 
that  the  null  hypothesis  of  no  cointegration  was  rejected.  Along 
with  the  Pesaran  critical  value,  this  study  was  also  used  the 
Narayan  [72]  critical  value  to  compare  the  calculated  F-statistics 
obtained  from  the  Wald  joint  test  of  significant  of  the  respective 
lagged  variables.  The  Narayan  [72]  critical  values  were  developed 
through  the  application  of  stochastic  simulations  specific  to  the 
sample  size  based  on  40,000  replications. 

However,  the  null  hypotheses  of  no  cointegration  were  still 
rejected  because  the  calculated  F  value  of  the  joint  test  (5.444)  was 
greater  than  the  upper  bound  critical  value  of  Narayan  [72]  (4.223) 
at  the  5%  significant  level  when  C02  emission  is  normalised.  Thus, 
the  long-run  cointegrating  relationship  among  the  respective 
variables  was  recognised  when  C02  emissions  were  normalised 
as  regressive.  Nevertheless,  when  the  GDP  growth  and  energy 
consumption  per  capita  were  considered  dependent  variables,  the 
calculated  F-statistics  fell  below  the  lower  bound  of  the  critical 
values  of  Pesaran  [69]  and  Narayan  [72], which  implied  that  no 
long-run  cointegration  existed  when  these  two  variable  were 
normalised.  Conversely,  when  population  growth  was  considered 
as  a  dependent  variable,  the  calculated  F-statistics  fell  between  the 
upper  and  lower  bounds  of  the  critical  value;  therefore,  the 
evidence  of  this  relationship  is  inconclusive. 

After  identifying  the  co-integration  relationship,  this  study 
preceded  to  estimate  Eq.  (3)  by  following  the  ARDL  specifications 
(1,  0,  0,  0,  1)  to  reveal  the  long-run  elasticity  of  the  respective 
variables  on  the  C02  emissions. 

The  estimated  long-run  coefficient  of  LGDPC  was  negative  and 
statistically  insignificant,  which  implied  that  the  C02  emissions 
initially  fell  with  a  rise  of  GDP  per  capita  as  shown  in  Table  5. 
However,  the  coefficient  of  the  quadratic  form  LGDPC  (LGDPCS) 
was  positive  and  statistically  significant,  which  indicated  that  the 
relationship  between  C02  emissions  and  economic  growth  was 
non-monotonic.  These  relationships  are  clearly  demonstrated  in 
Figs.  2  and  3,  in  which  C02  emissions  decreased  from  1970  to  1978 
and  then  smoothly  increased  from  1978  to  2012.  According  to 
Brock  and  Taylor  [73],  the  relationship  between  economic  growth 
and  environmental  quality  largely  depends  on  three  important 
mechanisms  such  as  scale  of  production,  composition  or  means  of 
production  and  use  of  technology  for  production.  In  relation  to  the 
study  findings,  during  the  period  of  1970-1978,  composition  of 
production  factor  was  less  responsive  towards  C02  emission, 
indicated  that  an  increase  of  GDP  per  capita  might  be  helpful  to 
reduce  per  capita  C02  emission.  However,  over  the  period  of  1978- 
2009,  the  positive  coefficient  of  LGDPCS  shows  higher  growth  rate 
of  C02  emission  than  GDP.  This  might  be  happened  as  agricultural 
sector  contributed  more  than  25%  of  Malaysian  GDP  from  1970  to 
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Table  4 

Results  of  the  bounds  tests  of  Eq.  (1). 


Dep.  var. 

SIC  lag 

F-statistic 

Probability 

Outcome 

Fco2(C02ILGDPC,  LGDPCS,  POPG,  EU) 

3 

5.444*** 

0.002 

Cointegration 

Flgdpc  (LGDPCIC02, LGDPCS,  POPG,  EU) 

3 

1.760 

0.165 

No  cointegration 

EpopctPOPG!  C02,LGDPC,LGDPCS,  EU) 

3 

2.671 

0.051 

Inconclusive 

Feuc(EUCI  C02,LGDPC, LGDPCS,  POPG) 

3 

0.648 

0.666 

No  cointegration 

Critical  value,  Pesaran  et  al.  [69] 

1(0) 

KD 

Narayan  [72]  /(0) 

KD 

1%  Significance  level 

3.29 

4.37 

4.280 

5.840 

5%  Significance  level 

2.56 

3.49 

3.058 

4.223 

10%  Significance  level 

2.20 

3.09 

2.525 

3.560 

Note :  Statistically  significant  at  1%  (***). 


Table  5 

Estimated  long  run  coefficients  using  the  ARDL  approach  (1,  0,  0,  0, 1 )  based  on  the 
Schwarz  Bayesian  criteria.  The  dependent  variable  is  per  capita  C02  emission,  and 
38  observations  were  used  for  estimations  from  1972  to  2009. 


Regressor 

Coefficient 

Standard  error 

T-ratio  [Prob] 

LGDPC 

-0.021 

0.013 

-1.630  [0.113] 

LGDPCS 

0.002** 

0.001 

2.337  [0.026] 

POPG 

0.001 

0.097 

0.012  [0.990] 

EUC 

0.235*** 

0.510 

4.607  [0.000] 

C 

0.869*** 

0.287 

3.019  [0.005] 

Note:  Statistically  significant  at  1%  (***),  5%  (**)  and  10%  (*)  level. 


Fig.  2.  Inverted  U  shaped  relations  between  per  capita  emissions  and  GDP. 


50  55  60  65  70  75 

Ipoly  smoothing  grid 


95%  Cl  -  Ipoly  smooth:  C02 


Fig.  3.  Inverted  U  shaped  and  upwards  trends. 


1978  (WD1,  2014)  and  then  share  of  agriculture  to  the  GDP  was 
decreased  while  shares  of  manufacturing  and  service  sector  were 
risen.  Al-Mamun  et  al.  [74]  also  revealed  that  economic  transfor¬ 
mation  from  industrial  to  service  sector  leads  to  reduce  C02 
emission  in  middle  and  lower  middle-income  countries.  Despite 
in  OECD  and  high-income  countries,  technologically  advanced  and 
energy-led  service  sector  increases  C02  emission  level  [74],  Like¬ 
wise  increasing  service-led  economy  emits  more  C02  emission 


Table  6 

Error  correction  representation  of  the  ARDL  model  (1,  0,  0,  0, 1)  based  on  SIC.  The 
dependent  variable  is  AC02,  and  38  observations  were  used  for  the  estimations 
from  1972  to  2009. 


Regressor 

Coefficient 

Standard  error 

T-ratio  [Prob] 

ALGDPC 

-0.011 

0.006 

-1.900  [0.066] 

ALGDPCS 

0.001 

0.655 

2.386  [0.023] 

APOPG 

0.662 

0.053 

0.012  [0.990] 

AEUC 

0.630 

0.186 

3.375  [0.002] 

AC 

0.477 

0.207 

2.296  [0.028] 

ecm  ( - 1 ) 

-0.549 

0.147 

-3.723  [0.001] 

ecm  —  C02 +0.02 1832  x  LGDPC - 0.0028490  x  LGDPCS-0.0012056  x  POPG- 
—  0.2354E-3  x  EUC-  0.86958  x  C. 


compared  to  other  sectors  [73,75],  Moreover,  some  studies  also 
demonstrated  that  higher  national  income  does  not  necessarily 
mean  a  greater  amount  of  pollutant  emissions  [76-79], 

The  findings  of  positive  relationship  between  economic  growth 
and  C02  emissions  are  partially  in  line  with  other  Malaysian 
studies  of  Ang  [61]  and  Azlina  and  Mustapha  [63]  although  they 
did  not  follow  an  EKC  path.  In  contrast,  Saboori  et  al.  [56]  found  an 
existence  of  EKC  for  Malaysia.  The  differences  of  this  study 
findings  may  occur  due  to  the  two  additional  control  variables, 
such  as  energy  consumption  and  population  growth  and  longer 
time  period  of  1970-2009  whereas  Saboori  et  al.  [56]  considered 
only  the  GDP  per  capita  and  its  quadratic  form  as  explanatory 
variables  for  a  period  of  1980-2009.  Furthermore,  measuring  the 
nonlinearity  by  taking  the  quadratic  form  fulfilled  the  necessary 
condition  for  the  existence  of  a  U  or  inverted  U  shaped  relation¬ 
ship  but  did  not  provide  the  sufficient  conditions,  as  indicated  by 
Sasabuchi-Lind-Mehlum  (SLM)  U  test. 

The  model  also  found  a  positive  but  insignificant  coefficient  of 
population  growth.  This  indicates  that  population  growth  was  not 
the  main  responsible  to  impact  the  C02  emissions  in  Malaysia. 
However,  energy  consumption  per  capita  had  a  strong  and  positive 
impact  on  C02  emission,  as  a  huge  proportion  of  C02  emissions  in 
Malaysia  comes  from  energy  consumption  of  the  power,  industrial 
and  transport  sectors.  This  relationship  is  in  line  with  many  other 
studies  e.g.  [10,20,61,80,81],  However,  technological  progress 
reduces  C02  emission  intensity  by  increasing  energy  efficiency 
[13-17]  where  carbon  free  new-technology  (e.g.  wind,  nuclear)  is 
conducive  in  reducing  C02  emission  without  harming  the  eco¬ 
nomic  growth  [82], 

Alike  the  long-run  relationship,  the  short  run  LGDPC  led  to  an 
approximate  1.1%  decrease,  whereas  the  LGDPCS  led  to  0.1% 
increase  in  C02  emissions  (Table  6).  In  addition,  Table  6  reports 
that  population  growth  does  not  foster  emissions  in  the  short  run. 
However,  in  the  short  run  as  expected,  the  energy  consumption 
per  capita  had  a  positive  and  significant  impact  on  the  C02 
emissions.  The  coefficient  of  the  equilibrium  correction  mechan¬ 
ism  (ECM)  was  0.549,  which  was  significant  at  the  1%  significance 
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Table  7 

ARDL-VECM  model  diagnostic  tests. 


R2= 0.83,  Adjusted  R2  =  0.80 

A:  Serial  Correlation /2(1)=0.012[0.912] 

B:  Functional  Form  j2(l)  =  1.828(0.176] 


C:  Normality  f(2)= 1.079(0.583] 

D:  Heteroscedasticity  j2(l)= 0.448(0.503] 


A:  Lagrange  multiplier  test  of  residual  serial  correlation;  B:  Ramsey's  RESET  test  using  the  square  of  the  fitted  values;  C:  Based  on  a  test  of 
skewness  and  kurtosis  of  residuals;  D:  Based  on  the  regression  of  squared  residuals  of  the  squared  fitted  values. 


Plot  of  Cumulative  Sum  of  Recursive 


Residuals 


Plot  of  Cumulative  Sum  of  Squares 
of  Recursive  Residuals 


The  straight  lines  represent  critical  bounds  at  5%  significance  level  The  straight  lines  represent  critical  bounds  at  5%  significance  level 


Fig.  4.  Plot  of  CUSUM  and  CUSUMQ  for  coefficients  stability  of  ECM  model  (critical  bounds  at  5%  significance  level). 


level  and  implied  that  disequilibrium  in  the  short  run  was 
adjusted  by  54%  per  year  towards  the  long-run  equilibrium. 

The  model  diagnostics  for  the  value  of  R 2  and  adjusted  R2  were 
83%  and  80%,  respectively,  which  indicated  that  the  model  was 
well  fitted.  The  ARDL  test  also  passed  several  diagnostic  tests,  with 
no  serial  correlations,  functional  errors,  normality,  or  heterosce¬ 
dasticity  problems  of  the  models  estimated  in  this  study  (Table  7). 

The  cumulative  sum  (CUSUM)  and  cumulative  sum  of  squares 
(CUSUMQ)  plots  (Fig.  4)  from  a  recursive  estimation  of  the  model 
also  indicated  that  the  model  was  stable  because  the  residuals 
were  within  the  critical  bounds  of  the  5%  significance  level. 


6.  Robustness  analysis 

6.J.  Dynamic  OLS 

The  robustness  of  the  long  run  coefficients  obtained  from  the 
ARDL  estimator  were  evaluated  by  applying  an  alternative  single 
equation  estimator,  such  as  the  dynamic  OLS  (DOLS)  procedure. 
The  prime  benefit  of  the  DOLS  approach  was  that  it  considered  the 
presence  of  a  mixed  order  of  integration  of  the  respective  variables 
in  the  cointegrated  framework.  The  estimation  of  DOLS  involved 
regressing  one  of  the  /( 1 )  variables  against  other  variables,  some  of 
which  are  /(l)  with  leads  (p)  and  lags  (-p)  of  the  first  difference 
and  others  are  1(0)  variables  that  include  a  constant  term  [75], 
Thus,  this  estimator  solved  two  important  limitations:  a  possible 
endogeneity  problem  and  small  sample  bias.  Moreover,  the  obtained 
cointegrating  vectors  from  DOLS  estimators  were  asymptotically 
efficient. 

The  result  from  DOLS  was  consistent  with  the  ARDL  according 
to  the  sign  and  significance  of  the  coefficient  as  presented  in 
Table  8.  The  negative  and  statistically  significant  coefficient  of 
LGDPC  indicated  that  C02  emissions  decreased  with  a  rise  of  GDP 
growth  during  the  first  phase  of  the  sample  years.  This  finding  was 
compatible  with  the  history  of  economic  growth  in  Malaysia.  From 
1970  to  1980,  the  primary  sector  mainly  agriculture  contributed 
42.7-33.0%  of  the  total  GDP  [6],  Based  on  this  fact,  the  primary 
sector  was  shown  to  be  much  less  responsible  for  the  C02 
emissions  than  the  manufacturing  sector,  which  is  reflected  in 
the  results  of  this  study.  However,  the  coefficient  of  the  quadratic 
form  of  the  LGDPC  was  positive  and  statistically  significant,  which 
indicated  that  during  1980-2009,  an  increase  in  GDP  growth 
intensified  the  C02  emissions.  This  finding  might  be  explained 


Table  S 

Results  from  the  dynamic  OLS. 


C02  emissions 

Coefficient 

Robust  std.  error 

z 

P-value 

LGDPC 

-0.0126*** 

0.0045 

-2.8200 

0.0050 

LGDPCS 

0.0029*** 

0.0003 

8.7100 

0.0000 

LEUC 

0.0002*** 

0.0000 

14.8600 

0.0000 

POPG 

0.0066 

0.0279 

0.2400 

0.8140 

Constant 

0.8479 

0.0744 

11.38 

0.0000 

R2:  0.967 

Note:  Statistically  significant  at  1%  (***),  5%  (**)  and  10%  (*)  level. 


by  Malaysia’s  economic  growth  from  1980  onwards,  when  the 
manufacturing  sectors  expanded  sharply.  The  transition  of  eco¬ 
nomic  growth  from  primary  to  secondary  sector  caused  a  sig¬ 
nificant  amount  of  energy  demands;  consequently,  the  trend  of 
C02  emissions  followed  a  U-shaped  direction.  Fig.  2  also  shows 
that  the  C02  emissions  increased  sharply  from  1980.  The  coeffi¬ 
cient  of  energy  use  was  still  positive  and  significant,  which  was 
consistent  with  the  long-run  coefficient  of  the  ARDL  estimators. 
However,  the  impact  of  the  population  growth  was  still  insignif¬ 
icant  on  the  C02  emissions  for  the  Malaysian  economy. 


6.2.  SLM  U  test  -  sufficient  condition  for  a  quadratic  relationship 

As  previously  mentioned,  the  estimation  of  the  ARDL  model 
takes  the  quadratic  form  of  the  LGDPC  that  fulfilled  the  necessary 
conditions  for  the  existence  of  a  U-shaped  or  an  inverted  U-shaped 
relationship;  however,  according  to  SLM  U  test,  these  conditions 
are  not  sufficient.  To  fulfil  the  sufficient  conditions,  the  SLM  U  test 
noted  that  the  conventional  econometric  model  is  no  longer 
suitable  to  test  the  composite  null  hypothesis  that  is  decreasing 
at  the  left  side  of  interval  relationship  and  increasing  at  the  right 
side  of  the  interval,  or  vice-versa.  Thus,  further  confirmation  of  the 
existence  of  a  U-shape  requires  to  perform  Sasabuchi-Lind- 
Mehlum  U  test  [66],  To  accomplish  this  task,  the  following  model 
should  be  estimated: 

LC02t  =  aGDPGt+bGDPG2  +ZcC+et  (5) 

The  joint  hypothesis  test  should  then  be  conducted 
H0  :  (a  +  b 2GDPGmin  <  0)  U  (a  +  b2GDPGmm  >  0) 
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Table  9 

Sasabuchi-Lind-Mehlum  test  for  U-shaped 
relationships. 


C02  emissions 

Slope  at  LGDPCm,n 

-0.497*** 

(-2.363) 

Slope  at  LGDPCmax 

1.458*** 

(2.384) 

SLM  test  for  inverse  U  shape 

2.36 

P-value 

0.0119 

Note:  Statistically  significant  at  1%  (***),  5%  (**)  and 
10%  (*)  level. 


This  should  be  conducted  against  the  alternative  hypothesis 

Hi  :  ( a  +  fa2GDPGmin  >  0)  U  (a  +  b2GDPGmax  <  0) 

where  LGDPCmin  and  LGDPCmax  represent  the  maximum  and 
minimum  value  of  GDP  growth.  If  the  null  hypothesis  is  rejected, 
this  confirms  the  existence  of  the  U  shape. 

Table  9  clearly  demonstrates  that  the  lower  bound  slope  of  the 
LGDPC  is  negative  (-0.497)  and  the  upper  bound  slope  of  the 
LHFC  is  positive  (1.458).  Both  bounds  are  statistically  significant, 
which  means  that  the  null  hypothesis  of  no  U-shape  is  rejected 
against  the  alternative  hypothesis  of  the  U-shaped  relationship 
between  C02  emissions  and  the  LGDPC.  The  negative  coefficient  of 
the  lower  bound  and  positive  coefficient  of  the  upper  bound 
endorse  the  validity  of  the  results  obtained  using  the  ARDL  and 
DOLS  approaches. 


7.  Conclusion 

This  study  found  several  interesting  findings.  First,  there  is  an 
existence  of  a  long-run  co-integrating  relationship  among  C02 
emissions  and  GDP  growth,  energy  consumption  and  population 
growth  in  Malaysia.  Second,  the  theory  of  an  environmental 
Kuznets  curve  is  not  valid;  instead,  a  U-shaped  curve  (opposite 
of  EKC)  exists  by  the  all  approaches,  such  as,  ARDL  bounds  test, 
dynamic  OLS  and  SLM  U-test.  To  check  the  stability  of  the  model 
variables,  the  cumulative  sum  (CUSUM)  and  cumulative  sum  of 
squares  (CUSUMQ)  plots  were  used  where  the  residuals  were 
within  the  critical  bounds  of  the  5%  significance  level.  During  the 
period  of  1970-1980,  per  capita  C02  emissions  decreased  with 
increasing  per  capita  GDP  (economic  growth);  however  since 
1980-2009,  per  capita  C02  emissions  increased  sharply  with  a 
further  increase  of  per  capita  GDP.  Finally,  the  analysis  reveals  that 
both  per  capita  energy  consumption  and  per  capita  GDP  has  a  long 
term  positive  impacts  with  per  capita  carbon  emissions,  but 
population  growth  rate  has  no  significant  impacts  on  per  capita 
C02  emissions.  However,  the  study  suggests  that  in  the  long  run, 
economic  growth  may  have  an  adverse  effect  on  the  C02  emissions 
in  Malaysia. 

Overall,  Malaysia  is  facing  an  increasing  trend  of  per  capita  C02 
emission,  GDP  and  energy  consumption  while  the  population 
growth  rate  shows  decreasing.  Thus,  the  rise  of  per  capita  C02 
emission  brings  a  huge  concern  to  the  increasing  demand  of 
energy  consumption  and  rapid  acceleration  of  economic  growth. 
Fossil  energy  still  possesses  absolute  share  of  the  total  energy 
consumption  although  Malaysia  has  adopted  the  Five-fuel  Diver¬ 
sification  Strategy  energy  mix  since  1999,  whereby  the  five  main 
energy  sources  are  oil,  natural  gas,  coal,  hydro  and  renewable 
energy.  Furthermore,  the  Government  of  Malaysia  in  2011  passed 
the  Renewable  Energy  (RE)  Act,  and  launched  a  roadmap  to 
increase  the  contribution  of  RE  to  the  electricity  generation  mix 
from  less  than  1%  at  present  to  5.5%  by  2015  [83]  Under  the  Act, 


since  December  2011,  a  Feed-in  Tariff  (FiT)  mechanism  is  imple¬ 
mented  to  be  paid  by  the  contributions  of  additional  1%  of  total 
electricity  bills  into  RE  fund  if  any  consumer  used  more  than 
300  kWh  of  electricity  per  month.  Instead,  this  mechanism  may 
work  as  an  incentive  for  energy  savings  and  a  motivation  to 
consumers  (individual  or  industry)  to  offset  the  incremental 
electricity  cost  by  applying  renewable  energy  and  energy  effi¬ 
ciency  measures.  Thus,  significant  transformation  of  low  carbon 
technologies  such  as  renewable  energy  and  energy  efficiency  could 
contribute  to  reduce  the  emissions  and  sustain  the  long  run 
economic  growth. 
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